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SUMMARY 

An analysis of two X-ray absorption methods for determining 
sTna.1 1 air-density gradients is presented. One method utilizes a 
Geiger-Mueller counter for detection and the other uses photographic 
film for detection. The methods are a refinement of a previous 
similar application. With the resulting sensitivity in density 
measurement , the first method was successfully applied to the evalu- 
ation of a 6-percent density gradient of a flat -plate boundary layer 
of about 0.024- inch thickness in a two-dimensional-flow channel. 

In the other method, the divergence of the X-ray beam passing the 
absorber poses a problem. A solution to this problem is postulated. 

In addition, example calculations of the expected sensitivity in 
boundary- layer density measurements are given for X-ray parameters 
selected for optimum operation. 


INTRODUCTION 

The application of soft X-ray absorption as a technique for the 
determination of air densities has a potential advantage over the 
optical radiation methods, especially in determining air densities 
over a long path length, in that X-ray refraction errors are reduced 
by at least a factor of 10^. In addition, all radiation methods have 
the advantage that no probes have to be inserted into the stream when 
evaluating density gradients in a wind tunnel. Radiation methods, 
however, can only be employed when two-dimensional flow can be pre- 
supposed. In references 1 and 2, the possibilities of the X-ray 
absorption method are recognized and are applied to the investigation 
of densities in wind tunnels - (reference l) and shock waves in super- 
sonic flow (reference 2) . All investigators using this method con- 
cluded that the full potential of the X-ray absorption method can only 
be realized if the probing X-ray beam can be made substantially smaller 
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than that used in references 1 and 2 by providing an X-ray source of 
higher intensity and an X-ray detector of higher sensitivity. 

A study of two different X-ray detecting methods was made at the 
NACA Lewis laboratory and the results and the conclusions obtained are 
presented herein. The sensitivity of the air-density measurement can 
be increased fourfold, even though the probing X-ray beam would be 
decreased in size by a factor of 30 with reference tc the instrument 
previously described (references 1 and 2) by employing Geiger -Mueller 
counters instead of an ionization chamber to provide a detector with 
increased sensitivity and by utilizing a tungsten-target X-ray tube 
instead of one with a chromium target to increase the X-ray intensity. 
With a resulting sensitivity in air density of ±1.5 percent at atmos- 
pheric pressure and a probe width of 0.002 inch, this method can be 
applied to density determinations in a boundary layer of an airfoil. 

The density gradient at one station in the boundary layer of a flat 
plate at subsonic flow of Mach number of 0.55 was experimentally 
determined by X-ray absorption measurements and compared with a cal- 
culated laminar boundary layer to indicate the feasibility of the 
X-ray-absorption method for such studies. 

Employing photographic film as an X-ray detector permits the 
analysis of a complete boundary- layer station from a single exposure 
without loss in sensitivity in air-density measurement. This usage 
eliminates the machining of accurate and movable parts, which are 
required in the point -by-point method. In the case of X-ray detection 
by photographic film, each individual ray acts as a probing element. 

The resolution of the measurement, however, is* determined by the 
geometry of the setup and by the size of the densitometer pickup used 
to measure the film densities. 

The X-ray beam exhibits a conical geometry; this divergence of the 
beam should pose no problem when employing Geiger -Mueller counters with 
an X-ray beam receiving slit of 0.002 inch width or a scanning X-ray 
beam with photographic film as the detector. For the single- exposure 
photographic method, however, the measured X-ray intensity must be 
corrected for the divergence of the X-ray beam. Suggestions for 
applying such corrections are presented. 

Operating characteristics required when using X-ray equipment with 
Geiger -Mueller counters and photographic film for air-density measure- 
ments in a boundary layer of an airfoil at given flow conditions are 
calculated. 
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THEORY AND DESCRIPTION CF METHOD 


The intensity I of an X-ray team remaining after absorption is 
given by Lambert ' s law 


I_ = e -BpL 
*0 


( 1 ) 


where 


I Q initial X-ray intensity 
B total mass absorption coefficient 


p density of gas 

L path length through absorbing gas 

(For convenience., all symbols are defined in appendix A.) This law 
presupposes that the initial X-ray beam is directly incident on the 
absorbing medium and that the source consists of parallel rays. If, 
however, the initial X-ray source is a point source and a distance 
d away from the detector, equation (l) becomes 


where the constant C is a function of the absorption that the initial 
X-ray beam experiences before entering and after leaving the absorber 
and A is the cross-sectional area of the X-ray beam at the distance 
d, as shown in figure 1. 

The sensitivity in air density is calculated from equation (2) as 


I_ = CA e "- BpL 


( 2 ) 


Io 4 ltd. 2 


dp _ dl 


( 3 ) 


p IBLp 


and. the minimum detectable density change can be expressed as 


Ap = 



( 4 ) 


BL 


which converts to equation (3) for 
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In this equation, Al/l represents the accuracy in intensity measure- 
ment . 

From equation (2), it is seen that the number of quanta entering 
the detecting device depends not only on the density of the absorber 
but also on the length of the absorbing path, the geometrical arrange- 
ment of the equipment, the cross-sectional area of the X-ray probing 
beam, the initial X-ray intensity, and the mass absorption coefficient. 

A study of these variables shows that: 

1. The absorber length is invariant in any given problem. 

2. The geometrical arrangement of the equipment should be such 
that the loss in X-ray intensity between the X-ray source and the 
absorber is a minimum. In the case where the mechanical geometry 
requires that the X-ray tube be mounted a given distance from the 
section containing the absorbing medium, this distance is kept as small 
as feasible and the loss in intensity due to absorption is minimized by 
passing the X-ray beam through an evacuated tube. Especially thin 
windows are employed at the entrance and exit sides of the absorber 
section for least intensity losses. In special cases, demountable 
X-ray tubes may be employed to altogether eliminate window material. 

3. The cross-sectional area of the probing X-ray beam is chosen to 
be a min imum for highest resolution but sufficient in size to give the 
intensity required by the detector for desired accuracy. 

4. A high initial X-ray intensity is desired to obtain a maxi mum 
X-ray flux for high detecting sensitivity and short time measurements. 
The initial X-ray intensity increases directly with the atomic number 
of the target material and the applied current and is a function of the 
applied voltage, as shown in figure 2. The tungsten- tar get X-ray tube, 
because of its higher atomic number, generates higher intensities than 
the iron- target tube. The intensity values in figure 2 were measured 
through the beryllium window of the X-ray tubes and at a fixed distance; 
hence the values are not absolute but rather are representative of 
relative intensities. The maximum value of the intensity is limited by 
the target material available in commercial X-ray tubes, by the current 
limitations set by the X-ray equipment, and by the maximum voltage that 
can be applied without loss in the desired sensitivity in the 
air-density measurement, as will be explained later. 

5. The mass absorption coefficient is selected to permit suffi- 
cient flux to reach the detector so as to approach the desired sensi- 
tivity in air-density measurement. The total mass absorption coeffi- 
cient for monochromatic X-rays, in this case for the effective wave- 
lengths of the generated continuous spectra, equals 
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.4 


(5) 


where S is a term due to scattering, which can he neglected because 
it never exceeds 0.1 percent for the soft X-rays employed in measuring 
air-density gradients and Z is the atomic number of the absorber, 
which is air in this study. Thus equation (5) indicates that the mass 
absorption coefficient is inversely proportional to the third power of 
the voltage, where the voltage is inversely proportional to the effec- 
tive wavelength of the applied X-rays. After the proportionality 
constant C' was experimentally obtained from two determinations of 
the mass absorption coefficient at 3.2 and 4.0 kilovolts, approximate 
mass absorption coefficients were calculated for several voltages and 
are plotted in figure 3. 

For high sensitivity in air-density measurements, according to 
equation (3) , a high mass absorption coefficient is desired and can be 
achieved by employing low voltages. A high intensity is also required 
to give sufficient counting accuracy and can be obtained by applying 
high voltages. The voltage is therefore adjusted to give sufficient 
X-ray intensity and a mass absorption coefficient for highest sensi- 
tivity in air-density measurement. 


Because the sensitivity of the density measurements is of prime 
importance for the evaluation of density gradients, its limitation is 
analyzed for the two selected X-ray intensity detecting methods of 
employing the G-M (Geiger-Mueller) counter and, subsequently, the 
photographic film as receiving element. 


Counter resolving time. - For measuring the air densities by means 
of X-ray absorption, G-M counters were used with a scaling circuit, 
which allows simple handling and produces stable operating conditions. 
Not only are G-M counters more sensitive than ionization chambers, but 
also their sensitivity does not change with a variation in wavelength. 
Each quantum, independent of the X-ray energy, entering the counter 
produces a pulse that after suitable amplification will trigger the 
scaling circuit. If, however, the quanta are following each other at 
shorter intervals than the resolving time T of the counter permits, 
a counting error equal to R m T , where is the measured counting 

rate, is introduced. The true counting rate n can be obtained as 
follows : 


SENSITIVITY OF DENSITY DETERMINATION 


Geiger-Mueller Counter as Detector 
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n 


Em 


1-E T 


m 


( 6 ) 


The counter resolving time varies for each counter, hut was experimen- 
tally found to he 5.7 X 10~® minutes on the average. This counter 
resolving time introduces an error of about 0.6 percent for a counting 
rate of 1000 counts per minute and is less for lower counting rates. 


Counting sensitivity. - The X-ray quanta entering the G-M counter 
are randomly distributed. The counting of such randomly distributed 
pulses introduces an error due to the law of statistics . This relative 
probable error AN/N is given as : 

AH _ 4. 0 • 67 n/ N+ % (y 

N N-N b 


where N is the total number of counts and N-g is the number of back- 
ground counts, which are always present. The number of counts due to 
the quanta emanating from the X-ray source is N - N B , where N B can 
be neglected when N»N B . This condition exists for a short -duration 
high -intensity count, such as is desirable for measurements of dynamic 
wind-tunnel conditions . The background count N B for a counting period 
of 1 minute was found to be about 40 to 50. For a duration of measure- 
ment of 1, 2, and 3 minutes and a counting rate of 500 counts per minute 
the relative probable error according to equation (7) is ±3.5 percent, 
±2.5 percent, and ±2.0 percent. Hence, 500 counts per minute seem a 
required minimum intensity to give sufficient counting accuracy over 
short enough measuring periods. 


Experimental method and results . - The accuracy in counting and 
the sensitivity in air-density measurement were experimentally investi- 
gated. A Machlett X-ray tube, Type A-2 Diffraction Tube, with a 
tungsten target and two separate 0.5-millimeter-thick beryllium windows 
was employed. The absorbing medium, air, was contained in a test 
section of 4-inch path length with a cellophane window (0.001 in. thick) 
on each side . The X-ray beam passed through the first cellophane window 
the absorbing air, and the second cellophane window to the measuring 
G-M counter. The air density in the test section was changed by evacu- 
ation with the aid of a vacuum pump. The measured X-ray intensity 
varied over a ratio of 1 to 3 for about one-half atmosphere change in 
density for an applied X-ray tube voltage of 3.2 kilovolts and a 
resulting mass absorption coefficient of 90 square feet per pound. In 
order to eliminate the errors in the X-ray-intensity measurements 
resulting from fluctuations in the X-ray emission, a reference count 
with a second G-M counter attached directly to one window of 
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the X-ray tube was obtained. The ratio of the measuring and reference 
counts n/n e is proportional to the intensity ratio in equations (l) 
and (2) . The probable relative error of the ratio N/W-p is composed 
of the errors of both counts and can be expressed as 



The experimental counting accuracies of three reference ratios each 
obtained from ten measurements with the calculated values for the same 
number of counts are compared in the following table : 


Experimental Calculated 

(averaged from 10 points) (from equations (7) and (8)) 


% 

% 

N 

W 

a(n/n r ) 

an e 

AW 

a(n/w e ) 

W R 

“W 


N 

n/w e 

40 

3,530 

1,850 

0.525 

±2.02 

±1.15 

±1.61 

±1.98 

120 

13,450 

7,090 

.525 

±1.52 

± .58 

± .82 

±1.03 

200 

19,600 

10,260 

.524 

±1.54 

± .54 

± .69 

± .87 


From the preceding table, it seems that a counting accuracy of 
±1.5 percent is about the best obtainable with the commercially 
available components. This accuracy is more than sufficient for the 
intended wind-tunnel applications. 

With the described test setup and a counting accuracy of about 
±2.0 percent, a sensitivity in air-density measurement of better than 
±1.0 percent is obtained at atmospheric pressure (equation (3)). This 
sensitivity seemed suitable for application of the method to the density 
measurements of a boundary layer. Such an application will be described 
later in this report. 


Photographic Film as Detector 

Theory of film detector . - The characteristic curve of a photo- 
graphic material, emulsion density as a function of the logarithmic 
exposure, is a straight line over a considerable portion. The equation 
of this line can be expressed as 


D = r (log ^ + log t - log i) 


0 ) 
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where 

D film density 
7 slope of line 

I intensity impinging on film, in this case X-ray intensity 
t exposure time 

i inertia of photographing material, constant of emulsion for fixed 
developing procedure 

Substituting equation (2) in equation (9) gives 

D = -0.434 TBLp + K* 


and 


P = 


0.434 7BL 


+ K" 


( 10 ) 


where 


and 


i 0 c 

K' = r(l°g ~2~ + log 1 - log 1 


K" = 


K' 


0.434 yBL 


Both K' and E" are constant for a given experimental arrangement and 
exposure time. From equation (10), the sensitivity in air-density 
measurement is 


dp _ dD 

p ~ 0.434 TBLp 


(ID 


and a minimum detectable change in air density can be expressed as 

An = — (12) 

p 0.434 yBL 


where AD represents the accuracy in emulsion-density measurement. 
Equations (ll) and (12) are very similar to equations (3) and (4), only 


7fi\Z 


MCA IN 2406 


9 


the accuracy in counting intensity is substituted by the accuracy in 
the film-density measurement for any constant 7 of the photographic 
material . 

Experimental method and results . - In figure 4 is shown the char- 
acteristic curve of Kodak film Spectrum Analysis No. 1 for 4.0-kilovolt 
X-rays using Kodak Developer D-19 and a developing time of 5 minutes. 
Other photographic materials and developing procedures were investigated 
but were found less satisfactory. Figure 5 illustrates the experimental 
arrangement. The X-ray beam from the tungsten-target Machlett Type A-2 
Diffraction tube passes through the 5.1-inch-long test section, through 
one 0.001 -inch-thick cellophane window, and impinges on the photographic 
film. In order to obtain the characteristic curve in figure 4, the 
energy (iXt) incident on the photographic film was changed by varying 
the exposure time t while the X-ray intensity I remained, constant . 
The emulsion densities within the straight portion of the characteristic 
curve range from 0.7 to 2.3. The slope of this line is 7 = 1.66. For 
each given set of X-ray intensities determined by the parameters of the 
experiment and the measuring conditions, the exposure time is so 
selected as to keep the resulting emulsion densities within the range 
of the straight portion of the characteristic curve . 

The setup in figure 5 was used to check equation (10) experimen- 
tally. The test section was evacuated to different air densities. A 
constant intensity X-ray beam of 4.0 kilovolts and 17.5 milliamperes 
was passed through the test section at each of the known air densities 
and the film was exposed for a predetermined time . The film densities 
were measured on a photographic analyzer with a pickup 0.063 inch in 
diameter. Three sets of measurements were made and are plotted in 
figure 6. The exposure time for each set of data was so adjusted as 
to keep the emulsion densities within the range of the straight portion 
of the characteristic curve . In order to normalize the three sets of 
data, equation (10) was changed to 

p _ p = (13) 

p m p 0.434 7BL 

where p m was arbitrarily chosen to be the maximum density measured for 
each set of data and D m is the emulsion density for the air 
density p m . The dashed line in figure 6 is calculated from equa- 
tion (13) for an experimental value of B of 46.6 square feet per pound 
and L of 0.425 foot. Figure 7 shows the photographic exposures corre- 
sponding to the experimental points of run 2 in figure 6. 
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In order to calculate the sensitivity of the air-density measure- 
ment, it is necessary to know the accuracy in the emulsion -density 
measurement. For this purpose a number of measurements were made at 
constant air densities, at fixed X-ray conditions and identical exposure 
times (about 45 sec) . Three sets of data obtained at three air 
densities are shown in figure 8. The emulsion densities plotted on the 
graph were measured on the densitometer. The probable error in D was 
found to be ±0.006 and to be independent of the emulsion-density value. 
The probable error in emulsion density D was measured to be ±0.005 by 
a microphotometer with a slit pickup of 0.0004-inch width, which 
indicates approximately the same accuracy. The smaller size pickup of 
the microphotometer makes this instrument better suited to measurements 
of density gradients over the extremely short distances encountered in 
boundary layers and also where high resolution in density gradient is 
required. A probable error in D of ±0.006 results in a sensitivity 
in air-density measurement of ±0.6 percent at atmospheric density 
(equation (ll)), which corresponds to a minimum detectable difference 

in air density of about 4 x 10~ 4 pound per cubic foot. 

With the resulting sensitivity in air-density measurement of 
±0.6 percent, the photographic method seems applicable to density 
determinations in a boundary layer. This arrangement has the advantage 
that a complete boundary -layer station is obtained from a single 
exposure, which eliminates the movable parts required for scanning in 
the previous method. 

APPLICATION TO WIND TUNNEL BY USING GEIGEE -MUELLER 
C0UNTEBS AS DETECTOBS 

Photographs of the tunnel, X-ray equipment, and detectors are 
shown in figure 9. Figure 10 is a cross-sectional drawing of the 
tunnel. A G-M counter directly connected to one window of a Machlett 
tungsten-target X-ray tube, Type A-2 Diffraction, measures the refer- 
ence counts. The X-ray beam from the second window is passed through 
an evacuated chamber, a cellophane window, the tunnel section of 4-inch 
path length, and a second cellophane window to the measuring G-M 
counter. A receiving slit 0.002 inch in width and 0.25 inch in length 
is attached to the front of the detecting G-M counter in such a way that 
it, together with the G-M counter, can be rotated to be alined parallel 
to the airfoil. Because of the extremely small size of the receiving 
slit, it must be perfectly alined with regard to the model and the air 
stream. The airfoil is a flat plate 1 foot in length mounted on a 
micrometer arrangement that moves it back and forth in relation to the 
0.002 -inch probing X-ray beam. The X-ray probing beam is stationary 
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the leading edge of the plate. The tunnel had its maximum operational 
flow at a Mach number of 0.55 and was operated at near atmospheric 
density. An applied X-ray voltage of 3.2 kilovolts was chosen. 

The calibration curve is given in figure 11. When modified from 
equation (2) for the reference method, the equation of this curve 
becomes 


JL _ ke~ B P L (14) 


where k combines those intensity and absorption factors for both 
counts, which remain constant during each measurement. The calibration 
curve is obtained by evacuating the tunnel to various densities, which 
are determined from pressure-probe manometer readings recorded simul- 
taneously with the counting-rate-ratio measurements. Once this curve is 
established for a given X-ray setup, it is necessary only to lo.cate one 
point for each X-ray run because only k changes, which results in 
merely a parallel shift of the calibration line. The point establishing 
the parallel shift is best obtained from a counting- rate ratio corres- 
ponding to the free-stream density p 0 , which is determined from 
pres sure -probe measurements. 

In order to compensate for the interference of the X-ray beam with 
the presence of the airfoil. X-ray- intensity curves were taken for each 
run without air flowing through the tunnel. In these compensating 
curves, the X-ray intensity is measured as a function of the distance 
of the plate surface from the X-ray beam. In a properly alined instru- 
ment, the intensity is unaffected when the plate surface is 0.001 inch 
and more removed from the slit. The experimental points in figure 12 
are calculated from X-ray intensities that are corrected from compen- 
sating curves measured immediately prior to or after each tunnel run. 

The plotted points in figure 12 are obtained from X-ray-absorption 
data, whereas the solid line is calculated for a laminar boundary- layer 
curve using the Blasius function and neglecting heat transfer. The 
theoretical error in the experimental X-ray absorption density measure- 
ments was calculated to be about dbO.8 percent, which corresponds to a 
minimum detectable difference in air density of 6 X 10”^ pound per 
cubic foot. Each point was obtained from three 1-minute determinations. 
The experimental points appear to establish a common trend, indicated 
by the dashed line, with an accuracy better than the calculated accuracy, 
This common trend of the X-ray density runs indicates the potentialities 
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of the X-ray absorption method for boundary- layer density evaluations 
even in cases similar to the present one, where the flow conditions 
are extremely unfavorable for any density- gradient determinations. 

The one station boundary- layer curve of density p against the 
distance y from the surface of the model was obtained by reading 
the X-ray intensity and by subsequently calculating the density p 
for each measured distance y. In order to eliminate the tedious 
procedure of data reading and point-by-point calculation, a recorder 
computer was developed that can plot the changing air densities 
directly as a function of a linear scale, which in the case of a 
boundary- layer determination would be the y- scale (see appendix B) . 

Operating characteristics required when using X-ray equipment with 
Geiger-Mueller counters and photographic film for air-density measure- 
ments in a boundary layer of an airfoil at given flow conditions are 
calculated in appendix C. 


EFFECT OF DIVERGENT X-RAY BEAM ON 
BOUNDARY- LAYER EVALUATIONS 

An X-ray beam consists of bundles of divergent rays due to the 
mode of propagation of X-rays. If the density is constant throughout 
the space encompassed by the X-ray beam, the absorption will be prac- 
tically the same for divergent and parallel rays. If a density 
gradient exists throughout the space covered by a divergent beam, 
however, each ray will traverse layers of different density. This 
problem is encountered in the investigation of boundary layers by an 
X-ray absorption method. 

The Machlett, Type A-2 Diffraction, tube has a focal spot of 
about 0.004 by 0.04 inch at the target as seen from the direction of 
the beryllium window. In the X-ray arrangement with the G-M counter 
as detector, the 0.002- inch slit sees the 0.004- inch wide X-ray beam 
entering the test section with a width of less than 0.003 inch anr) as 
if the X-rays would emanate from an array of point sources. In 
addition, an integrated intensity over the receiving- slit area 
represents each point measurement . The divergence of the rays there- 
fore hardly affects the measured absorption in the case of point-by- 
point determinations provided the size of the detecting pickup is 
of the same order of magnitude as the beam probing the absorbing 
medium. 
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As for the single -exposure photographic method, the divergence of 
the team poses a problem due to the large area of the detector of more 
than 0.25 inch in diameter. In addition, an integrated intensity over 
the exposure is no solution because all density changes to be investi- 
gated are contained in this area. The individual measurements are 
obtained by integrating over small areas, which are determined by the 
size of the densitometer pickup. 


On the assumption that all rays diverge radially from a point 
X-ray source, which is directly adjacent to the entrance side of the 
test section, so that d = L, the air density through the traversed 
distance r a is integrated by each X-ray at a given angle cp and 

C p, C7 

0 from the center ray of the divergent beam. Then the emulsion density 
is a function of r, cp , and 9 and equation (10) for any constant 
cp and 9 becomes 

E 


angle 


Dcp,0 = - 0-434 7 f Bp(rcp^ 0 ) dr^e + K' 


(15) 


Because the experimental emulsion densities D^ q of the exposure are 
obtained as a function of the linear distance y from the model, it is 
desirable to obtain equation (15) in the same coordinates. For 
two-dimensional flow, P( r q 5 q) can substituted by p(y), 
r cp 0 = (y- a )/ sinC P 008 and cLqp q = dy/sin cp cos 0 , where a is 

the distance from the plate to the center of the X-ray beam, which is 
the center of the exposure. The length E of all r cp, 0 can assumed 

to be constant and equal to the path length L; then 


cos 0 sin cp 


Y - a 


L 


where Y is the distance of r^ q normal to the model at the location 
of the photographic material and equation (15) can be written 

D Y = ” ° - f- ' F + K ( 16 ) 

a 


Differentiation with respect to the parameter Y gives 


dD 

dY 


0.454 7 L 
(Y - a ) 2 


Bp(y)dy 


0.454 7 LBp(Y) 
Y - a 


(17) 


In order to determine the density p(Y), 
are combined to yield 


equations (16) and (17) 
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P(Y) 


1 

0.434 7BI 



x <*D 

a ) dY + D - K' 


(18) 


For the boundary condition, Y = a and p(Y) = p(a) . Then equation (18) 
reverts to the parallel -beam relation given in equation (10); 


therefore 


K' = D a + 0.434 7 BLp(a) (19) 

The photographic exposure gives an emulsion density against distance 
(D against Y) curve . The derivative of that curve dD/dY is required 
to obtain the air-density distribution through the boundary layer, 
(equation (18)) . The tangent cffi/dY can be measured for each point by 
mechanical means and equation (18) can thus be plotted point by point. 
The accuracy of this procedure depends entirely on the scale to which 
the emulsion -density curve D against Y is drawn; therefore the 
boundary -layer curve, equation (18), can be obtained as accurately by 
this means as the experimental data permit. An empirical curve -fitting 
procedure would appear to be another plausible method of attacking the 
problem. For instance, by determining a differentiable power relation 
for the experimental emulsion density and distance curve, equation (18) 
can be used to calculate point-by-point air-density values . 

When the X-ray source is not directly adjacent to the test section, 
so that d> L (fig. l), the density-distance equation becomes 

p(Y) = | p(Y) iQ + (1 - |) p(Y x ) (20) 

00 

- I p(Y) 18 + \ X (1 - I* 11 p(Y n>18 < 21 > 

n=l 

where is the density as given in equation (18) and p(Y^) is 

the density at the entrance points Y^ where Y^ is the distance 
measured at the X-ray entrance side to the absorber corresponding to Y. 
The entrance points Y^ are related to the measured distances Y so 

that 


Tl « (1 " j)Y + £ a 
d d 

The term p(Y n ) denotes the densities as calculated from equation (18) 
for points T l> Y 2* • • * V where 
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Y n - l 1 - j)Vl + (l) a 

Because the second factor of equation (21) in most cases is a rapidly 
converging series, the calculation of two to three terms will usually 
suffice. 


CONCLUDING REMARKS 

Commercial X-ray equipment with Geiger-Mueller counters as 
detectors was found adequate for measuring the densities within an 
accuracy of 6 X 10"^ pound per cubic foot in a tunnel of 4 inch width 
operating at near atmospheric density, where an X-ray-beam width of 
0.002 inch and a measuring time of about 3 minutes were used. For a 
1-minute measuring time, therefore, an accuracy in density measurement 
of 1.1 X 10- 3 pound per cubic foot was obtained. 

The X-ray technique of measuring density has been applied to a 

3 

boundary layer of a flat plate at a station 7— inches from the leading 
edge and at a Mach number of 0.55. 

For lower-density tunnels, higher X-ray intensities and softer 
X-rays are required to maintain the accuracy of the measurements. 

These requirements can be achieved by employing specially built, 
demountable X-ray tubes and by eliminating all window material. 

For more efficient operation and production of a permanent record 
of the density measurements, a counting rate computer recorder was 
developed. Experiments have shown this instrument to be successful in 
plotting air densities directly as a function of a linear scale. In 
order to record the air densities through a boundary- layer station, the 
motion of the paper has to be coupled to the motion of the probing 
X-ray beam. 

Commercial X-ray equipment with photographic film as a detector 
was found adequate for measuring the density in a 5.1-inch path length 
with an accuracy of about 4 X 10" ^ pound per cubic foot when a 
densitometer pickup of 0.063-inch diameter was used, and with a 
microphotometer slit of 0.0004 -inch width. An exposure time of less 
than 1 minute was employed in both cases. 

Photographic material as the detector has the advantage over 
counting methods in that all desired data are obtained from one meas- 
urement. In addition, no movable parts are required for scanning. 
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The photographic method poses a problem, however, in that a 
correction for the divergence of the X-ray beam seems necessary when 
evaluating density gradients such as those encountered in boundary 
layers. Correction methods are suggested. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, March 7, 1951. 
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APPENDIX A 
SXMBOLS 

The following symbols are used in this report: 

A X-ray beam area at distance d 

a distance from airfoil to center of X-ray beam 

B total mass absorption coefficient 

C function of constant absorption in d - L 

C ' proportionality constant 

T film emulsion density 

D film emulsion density at a 

a 

film emulsion density for p = P m 
film emulsion density at Y 

film emulsion density at X-ray angle cp and 9 
distance from X-ray source to detector 
X-ray intensity after absorption 

X-ray intensity at n kilovolts, n = 3.2, 4.0, 4.9, 6.8 
I Q initial X-ray intensity 

i inertia of film emulsion 

K proportionality constant, equation (5) 

K 1 proportionality constant, equation (10) 

K" proportionality constant, equation (10) 

k proportionality constant, equation (14) 

L absorber path length 

M Mach number 




D c 

d 

I 

I 


n 


j 
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N 

total measuring counts 

n b 

total background counts 

W E 

total reference counts M 

3>! 

n 

ro 

true counting rate 

E 

path length of diverging X-ray, E » L, E * d 

Em 

measured counting rate 

r,cp,<9 

spherical coordinates with X-ray source at center 

r cp,0 

path of X-ray divergent at angle cp and 0 

S 

absorption due to X-ray scattering 

t 

film exposure time 

V 

X-ray accelerating voltage 

Y 

distance from airfoil as measured on photographic film 

Y 1 

distance measured at X-ray entrance side to absorber corre- 
sponding to Y 

y 

distance from airfoil 

z 

atomic number 

7 

slope of characteristic film curve 

P 

air density 

p o 

free -stream air density 

P B 

air density at surface of model 

p m 

maximum air density for each run in figure 5 

P(Y) 

air density as function of Y 

T 

resolving time of Geiger -Mueller counter 
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APPENDIX B 
EECOEDEE COMFUTEE 

The X-ray intensities are usually obtained as counts per unit 
time . Commercial counting -rate meters that record the integrated 
pulses over a fixed time are available; most commercial counting -rate 
meters , however, are unsuitable for problems where rapidly changing 
counting rates of randomly distributed pulses have to be recorded, 
because a relatively long time determination at each constant condition 
is necessary if the measurement is to be accurate within the relative 
probable error as given in equation (7). An instrument was therefore 
developed that measures the time required for a constant number of 
counts instead of the counting rate . This instrument has the added 
advantage that the relative probable error introduced by the statis- 
tical nature of the occurrence of the pulses remains fixed provided 
the background count is small compared with the measured count. The 
air density is then a direct function of the logarithm of the time 
required to permit the passage of the fixed number of pulses. For the 
developed apparatus, the amount of time for a given air density is 
proportional to the angle of rotation of a cylinder that carries a 
logarithmic spiral. The cylinder is stopped in its rotation almost 
instantaneously at the occurrence of the last of the predetermined 
number of pulses. Provisions are made for printing a record where 
the spiral touches the paper chart. After printing, the cylinder, 
returns automatically and rapidly (in a few seconds) to its starting 
position to be ready for the next determination. 


Air densities measured by the X-ray absorption method have thus 
been recorded and check well with the air densities calculated from 
pressure and temperature readings. When recorded X-ray densities are 
calculated back to intensities, they are practically in complete 
agreement with counting-rate readings taken simultaneously with the 
recordings. Hence, the recorder does not seem to introduce any 
noticeable additional inaccuracies into the air-density measurements . 

For the experiments the paper chart moved on a ratchet arrangement . 
For boundary -layer evaluations, the motion of the paper will have to be 
coupled to the X-ray scanning mechanism. 
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APPENDIX C 
SAMPLE CALCULATION 

The practicability of the use of the X-ray-absorption method for 
air-density-gradient studies is demonstrated by the presentation of a 
typical boundary-layer problem encountered at the Lewis laboratory and 
the manner in which it would be treated using in turn the Geiger- 
Mueller counter and photographic film as receiver. 

A supersonic wind tunnel embodying a flat-plate model has the 
following operating parameters: 


Tunnel width, L, ft 3 

Free-stream density, pQ, lb/cu ft 0.0374 

Mach number, M 2.38 

which gives a boundary-layer density ratio, P b /Pq 0.47 

yielding a value for the density at the model 

surface, p B , lb/cu ft 0.0176 


Employment of the G-M counter necessitates calculation of 
the minimum applied X-ray accelerating potential to yield a sufficient 
counting rate, at least above 500 counts per minute, at the maximum 
density of the boundary layer (Pq) • 

A preliminary calculation using equation (2) and previously 
obtained data revealed that a mass absorption coefficient of approxi- 
mately 30 square feet per pound is needed for the desired counting rate 
at the given operating conditions of the 3-foot tunnel. A study of 
both figures 2 and 3 indicates that a voltage of 4.9 kilovolts, in a 
stable operating region, corresponds to a mass absorption coefficient 
of 27.0 square feet per pound. Hence a sample calculation is shown for 
an accelerating potential of 4.9 kilovolts. 

Under previously described experimental conditions and for an 
X-ray voltage of 3.2 kilovolts, 

I 3.2 - — ~2 CAe = 1700 (eounts/min) 


for 

A = ( 0. 002) x (0.25) (sq in.) 

d =» 9.1 (in.) 

B = 90 (sq ft/lb) 
p = 0.0730 (lb/cu ft) 

L - 4.0 (in.) 
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The relative intensity ratio I 4 g/lj 2 determined from fig- 
ure 2 is about 9. Then under the same conditions, 


E = -±°- 
l 4 . 9 . ,2 

4jtd^ 


CAe“® p -* J = 9X Ij g = (9)x(l700) (counts/min) 


yielding 


( x o H ) 


9 X 1700 d' 


2 9 X 1700 


0 / 4 . 9 -BpL 




0.112 


= 78,000 (counts/min) 


for the specified tunnel conditions of 

d = 41.1 (in.) 

B = 27.0 (sq ft /lb) 

L = 3.0 (ft) 

equation (2) at 4.9 kilovolts becomes 

= 18,000 e -81p = 6660 e -81p 


4.9 




For the maximum density, p^ = 0.0374 (lb/cu ft) 


I 4 9 = 6660 e -3 * 03 = (6660)X (0.049) 

= 326 (counts/min) 

Because this number of counts is below that minimum previously spec- 
ified, the use of a larger slit is suggested. A l/2-inch slit gives 
about 650 counts per minute and is assumed for the following calcu- 
lations . 
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For the minimum density p B = 0.0176 pound per cubic foot, 

I 4#9 = (2) 6660 e"°* 143 = (2)x(6660)x(0.867) 

= 11,560 (counts/min) 

For the lowest counting rate of 650 counts per minute and a ref- 
erence counting rate of 5000 counts per minute, the statistical error, 
calculated from equations (7) and (8), is approximately ±3.1 percent) 
hence 



±0.031 


For the given boundary layer from p^ = 0.0374 pound per cubic foot 
to pg = 0.0176 pound per cubic foot using a 0.002- by 0.5-inch receiv- 
ing slit, the minimum sensitivity is calculated from equation (3). 



±0.031 

27(0.0176)X(3) 


= ±2.2 (percent) 


The density change over the entire boundary layer (p Q -pg) is 

0.0198 pound per cubic foot) then the detectable density change, as 
calculated from equation (4), is less than ±2 percent of the entire 
density range of the boundary layer. 


In order to determine the feasibility of using the X-ray photo- 
graphic method for the given tunnel, the optimum combination of sen- 
sitivity and exposure time must be calculated. The mass absorption 
coefficient B best suited for the given operating conditions can be 
obtained from equation (12), where in this case ZD is the useful 
range in emulsion densities from 0.7 to 1.1 if the microphotometer is 
to be used for the measurements of emulsion density. 


Then 

B = (p - p ) 0 . 434yL = (0.0198)°(0.434) r L = 9,3 ^ S<1 ft / lb ) 

13 

Examination of figures 2 and 3 shows that this mass absorption 
coefficient corresponds to an accelerating voltage of about 6.8 kilo 
volts and is in a stable operating range for the X-ray equipment. 
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The exposure 

the energy (ixt) 
density. Thus 


time tg.g is calculated by utilizing the fact that 
remains constant for a fixed level of emulsion 


t 6»8 - H.O (T ) 


^oU.o 


0 y 6.8 


And making use of equation (2), 


, (l oh.O ,-(BpL) 

%.8 - *4.0 (I 0 ) 6 . 8 WoJ -<BpL) 


4.0 

6.8 


where the subscripts 4.0 and 6.8 refer to the parameters used with the 
experimental test chamber and the given tunnel, respectively. The 
initial X-ray-beam intensity ratio (l Q ) 4 q /^ q^q 8 = which is 

obtained from figure 2. The exposure time t 4< Q has been experimen- 
tally determined, as shown in figure 4, to be 30 seconds for D = 0.7. 

A substitution of the given numerical values into the preceding 
equation leads to the following result for the exposure time: 

2 -(46.6) (0.073) (-ri) 

t =30 (J £ 

6.8 \7 • 5 / \9 . 1 / -(9.3) (0.0374) (3) 

e 


= 55 (sec) 

The minimum sensitivity is determined from equation 

i 

dp\ dD _ ±0.006 

v p ) . ~ 0.434 y-BLp (0.434) (1.66) (9.3) (0.0176) (3) 
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= ±2.1 (percent) 


2 . 


24 


NACA TN 2406 



Figure 1. - Schematic diagram of setup. 
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Figure 2. - Variation of X-ray-beam intensity with applied tube 
voltage. Plate current, 17.5 milliamperes . 
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Figure 3. - Variation of mass absorption coefficient of air with voltage applied to 

X-ray tube. 
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Figure 4. - Characteristic curve for Kodak film Spectrum Analysis No. 1. X-ray voltage, 
4.0 kilovolts; current, 17.5 milliamperes; air density, p, 0.073 pound per cubic foot. 
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Figure 5. - Schematic diagram of X-ray photographic arrangement. 
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Figure 6. - Difference in air density as a function of emulsion 
density. Probable mean error in D - Dm, ±0.008, 
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Figure 7. - Air density against emulsion density. (Representation of run 2 from fig. 6.) 
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Figure 8. - Emulsion-density reproducibility at 4.0 kilovolts, 17 . 5-milliampere plate current, and 

45 -second exposure. Average AD, ±0.0055. 
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(a) X-ray beam entrance side. 

Figure 9. - Tunnel section with X-ray equipment and associated parts. 
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(b) X-ray beam exit and receiving side. 

Figure 9. - Concluded. Tunnel section with X-ray equipment and associated parts. 
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Figure 10. 


Cross section of tunnel. 
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Figure 11. - Calibration curves at X-ray voltage of 3.2 kilovolts and 

tunnel width of 4 inches. 
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Figure 12. 


Variation of density ratio with distance from model. 

number, 2.46 x 10 6 . 
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